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摘 要 本 文 以 内 蒸 蒸 发 机 理 的 特殊 性 为 出 发 点 ,基于 过 热 液 滴 闪 燕 所 需 热量 由 其 内 部 过 热能 量 提 供 这 一 理论 基础 ， 建 
立 了 全 新 的 过 热 水 滴 内 蒸 的 现象 学 模型 , 使 用 MATLAB/Simulink 进行 求解 ,通过 实验 验证 其 有 效 性 . 经 计算 得 出 了 内 


燕 过 程 中 液 滴 温 度 、 直 径 、 质 量 、 营 发 速率 和 蒸发 率 的 变化 规 得 


E， 并 利用 控制 变量 法 ,研究 初始 温度 、 闪 燕 压力 和 液 滴 粒 


径 对 液 滴 温度 变化 、 蒸 发 速率 、 蒸 发 率 和 内 燕 时 间 的 影响 ,得 出 了 具有 代表 性 的 结论 . 
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Abstract Considering the particularity of i of flash evaporation, a new phenomenolog- 
ion was built. This model was based on the idea that 


ical model for superheated droplet flash e 
the heat needed for evaporation was ided by the superheat energy inside the droplet. The model 


was calculated using MATLAB/Si 


ink and was validated by experimental data. 


The evolution 


of the droplet’s 机 地 iiaeter，Inass， evaporation rate and flash efficiency was obtained. 


The influence of initial tt Wture, evaporation pressure and droplet diameter on the droplet’s 
temperature evoiotion dation rate, flash efficiency and flash duration time was investigated and 
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0 引言 


液 滴 和 喷雾 的 蒸发 因 在 工程 上 有 重要 应 用 价值 
而 成 为 数 干 年 来 的 研究 热点 。 其 中 ， 喷 雾 内 燕 作 ， 
一 种 特殊 的 营 发 形式 , 在 工业 余热 利用 、 航 空 航 入“ 
内 燃 机 等 领域 的 应 用 前 景 广阔 ， 近 年 来 逐渐 受到 重 
视 . 液 滴 闪 燕 是 喷雾 内 燕 过 程 的 重要 组 成 部 分 ， 其 
燕 发 特性 和 传 热 传 质 机 理 具 有 较 高 的 研究 价值 

液 滴 内 燕 的 实验 和 理论 研究 均 取 得 一 定 进 
展 . Owen 等 日 对 直径 为 1~3 mm 的 液 滴 在 不 同 
过 热度 下 的 内 营 强 度 进 行 了 实验 研究 ， 发现 过 热 
度 对 液 滴 闪 花 强 度 有 重要 影响 . 周 致 富 等 四 以 制 
冷 剂 R134a 闪 蒸 喷雾 为 研究 对 象 ， 得 出 了 液 滴 直 
径 和 液 滴 速 度 在 喷雾 场 中 的 分 布 。Shin 等 加 首 
先 提出 了 扩散 控制 蒸发 模型 ， 用 于 预测 液 滴 真空 
闪 蒜 结 冰 过 程 ， 该 模型 在 液 滴 真空 办 蒸 研究 中 得 
到 广泛 应 用 3。 Cheng 等 9 在 此 基础 上 提出 了 
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\ 效 热 导 率 模型 ， 得 出 了 比 Shin 模型 更 为 准确 的 
。' 结 


论 。Muthunayagam 等 中 采用 蒸汽 扩散 模型 预 
测 液 滴 真空 内 莱特 性 ， 并 以 实验 验证 模型 的 有 效 
性 . Adachi 等 四 提出 了 过 热 液 滴 闪 蒸 过 程 中 传 热 系 
数 的 关联 式 ， 在 很 大 的 过 热度 范围 内 均 可 使 用 ， 因 
此 被 普遍 应 用 于 燃料 喷雾 闪 蒸 研究 领域 810. Raju 
等 HH 对 一 个 过 热 燕 发 模型 进行 改进 , 使 之 同时 适用 
于 过 热 和 非 过 热 蒸发 过 程 ， 并 使 用 Monte Carlo 方 
法 对 喷雾 办 蒸 进行 计算 ， 

综合 来 看 ， 前 人 对 于 液 滴 闪 蒸 的 研究 多 集中 于 
压力 极 低 、 接 近 真 空 的 范围 ， 所 用 工 质 多 为 液态 燃 
料 ; 所 选用 的 模型 ， 例 如 前 文 提 到 的 扩散 控制 蒸发 
模型 、 有 效 热 导 率 模 型 和 蒸汽 扩散 模型 等 ， 均 是 基 
于 蒸发 过 程 由 扩散 作用 驱动 这 一 假设 .而 在 工程 实 
际 中 ,过 热 的 水 在 常 压 甚至 高 于 大 气压 力 下 的 闪 蒸 
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亦 有 重要 的 应 用 价值 这 一 过 程 由 传 热 驱动 , 前 述 。 dla 4os(- 区 ) -6a (7 (3) 
模型 并 不 适用 . 因此 ,本文 运用 质量 守恒 、 能 量 守恒 。 于 2 Fh 
和 传 热 传 质 等 基本 定律 ， 建立 过 热 液 滴 在 较 高 压力 。 式 中 : cy 是 液体 比 定 压 热 容 ， pa 是 液体 密度 ，4 
下 闪 雍 的 现象 学 模型 , 研究 其 蒸发 过 程 和 影响 因素 ，” 液 滴 表面 积 , V 是 液 滴 体 积 
拓展 液 滴 内 蒸 的 研究 范围 ， 为 过 热 液 滴 内 蒸 的 工程 液 泣 质 量 的 变化 是 ， 
应 用 提供 参考 ， 
Var hv FE Ee ,EH 一 ~ 一 mov (4) 

1 过 热 液 滴 办 蒸 数学 模型 . 

本 文 所 研究 的 过 热 液 消 闪 燕 的 基本 过 程 是 , 初 。 式 中 ma 是 液 消 质 量 。 
温度 为 Tio、 初 始 直径 为 Dao 的 过 热 水 滴 ,在 压力 RE 
为 pev 的 水 兹 汽 环 境 中 发 生 办 燕 ， 环 境 温度 维持 在 Da- (om 加 
该 压力 对 应 的 饱和 温度 7。 为 了 简化 计算 过 程 ， 突 pa 

轩 蔡 也 考点， 设 : 1) 液 滴 在 闪 菊 六 Po 
出 内 敬 的 基本 特点 ， 作 如 下 假设 : 1) 液 滴 在 闪 敬 过 0 


程 中 保持 为 球形 ; 2) 将 液 滴 分 为 内 部 和 表面 两 个 区 ， 
认为 内 部 温度 均匀 变化 ， 没 有 温度 梯度 ;内 蒸发 生 
在 液 滴 表面 ， 表面 温 度 维 持 在 沸点 ，3) 忽略 液 滴 与 
周围 环境 的 辐射 及 对 流 换 热 ; 4) 不 考虑 液 滴 内 部 因 
过 热 产 生 的 微 爆 效应 ; 5) 忽略 液 滴 内 部 气泡 对 藻 发 


恒定 
闪 燕 与 普通 蒸发 在 机 理 上 最 主要 的 
通 燕 发 所 需 热量 主要 由 来 自 周围 气体 
字 ， 而 关 闪 大 主 要 由 请 能 坟 
供 , 燕 发 速率 受 液 滴 由 内 向 外 的 咎 秘 传 递 控制 . 在 较 
低 的 过 热度 下 ， 液 滴 内 部 的 时 头 和 对 流 是 主导 的 伟 
热 过 程 ; 在 很 高 的 过 热度 下 , 莱 汽 空 化 占 主 导 . 相关 
过 程 极其 复杂 , 无 法 建立 分 析 解 .因此 ,在 建立 液 滴 
闪 蒜 数学 模型 的 过 程 中 , 使 用 基于 经 验 的 传 热 系数 ， 
忽略 液 滴 内 部 气泡 形成 、 生 长 和 破碎 等 因素 的 影响， 
在 满足 工程 应 用 精度 的 前 提 下 ， 达 到 简化 计算 、 提 
高 效率 的 目的 . 闪 葵 燕 发 速率 表达 式 是 ， 
aa 和 


-Th _ pao 
h 

式 中 : mov 是 燕 发 速率 , ra 是 液 滴 半 径 ，7a 是 液 滴 

温度 , Tb 是 沸点 温度 ,hts 是 汽化 洪 热 , as 是 从 液 滴 

内 部 向 表面 的 总 传 热 系数 , 单位 是 kJ.s-1m?.K-1。 

该 传 热 系数 取 自 Adachi 提出 的 关联 式 回 ， 考 虑 了 

液 滴 内 部 成 核对 换 热 的 强化 ,涵盖 了 很 大 的 过 热度 


Tov = 4nr3as 
fg 


范围 : 
0.76AT026 (0< AT<5) 
as 一 4 0.027AT?33 (5<AT<25) (2) 
13.8AT0.39 (AT > 25) 


式 中 AT 是 液 滴 的 过 热度 , AT = Tu 一。 
液 滴 的 温度 变化 可 以 写作 : 


ET NA) 


论 - (mao = ma) /mao x 100% (6) 
式 站 | 了 是 液 消 的 蒸发 素 ， mao 是 深沉 的 初始 质量 - 
NS 求解 方法 
使 用 MATLAB 软件 的 仿真 平台 Simulink, 对 过 
热 液 滴 办 疗 过 程 进 行 建 模 与 求解 . 根据 方程 (1)~(6)， 
建立 如 图 1 所 示 的 系统 模型 为 使 系统 清晰 简洁 ， 
过 热度 的 指数 项 和 物性 的 表达 式 均 封 装 在 子 系统 内 . 
采用 固定 步 长 三 阶 Runge-Kutta 算法 , 通过 计算 得 
出 闪 燕 过 程 中 液 泣 温 度 、 吉 径 、 质 量 、 燕 发 速率 和 莱 
发 率 的 变化 情况 (人 
3 0 
有 效 性 验证 
人 人 验证 妆 学 模型 的 有 效 性 ， 在 课题 组 措 建 的 喷 
CN 台 上 进行 实验 ， 以 工程 中 重点 关注 的 参 
燕 发 率 为 研究 对 象 , 将 实验 结果 与 计算 结果 
es a 
测 ， 但 考虑 到 喷雾 闪 葵 的 基础 是 液 泣 闪 燕 ， 其 总 体 
的 蒸发 率 能 够 代表 单个 液 滴 办 燕 的 结果 ， 因 此 这 种 
验证 方式 是 可 行 的 .实验 系统 如 图 2 所 示 ， 该 实验 
台 的 特点 是 可 以 在 较 高 的 液体 温度 和 闪 薰 压力 下 安 
全 稳定 运行 ， 革 发 率 通过 计算 闪 节 蒸汽 流量 与 供水 
流量 之 比 获 得 . 
选取 3 组 不 同 的 工 况 ,最终 莱 发 率 的 计算 结果 
与 实验 结果 的 对 比 见 表 1. 可 以 看 出 , 计算 结果 与 实 
验 结果 相差 非常 小 ， 且 计算 结果 均 偏 低 ， 偏 差 的 可 
能 成 因 包 括 数学 模型 选用 的 经 验 公式 的 准确 度 、 计 
算 过 程 中 物性 的 准确 度 等 这样 的 误差 属于 工程 上 
可 以 接受 的 范围 ， 可 以 认为 本 文 建立 的 数学 模型 能 
够 准确 的 描述 液 泣 闪 蒸 过 程 ， 
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Fig. 1 2 model 


算得 出 的 上 述 各 参数 变化 曲线 如 图 3~7 所 示 . 可 以 
看 出 , 菊 发 在 刚 一 开始 速度 非常 快 , 液 滴 的 温度 、 直 
径 和 质量 都 迅速 降低 ， 其 中 液 滴 温 度 在 0.44 s 内 就 
下 降 了 90%。 随 着 液 滴 温度 趋 近 饱 和 温度 ， 藻 发 速 
率 迅 速 减 小 .与 之 相对 应 ， 液 滴 营 发 率 在 开始 后 较 
短 时 间 内 即 太 幅 中 着 至 接近 最 大 值 ， 此 后 缓慢 增加 
直至 不 再 改 礁 ;总体 来 看 ， 液 消 闪 六 在 很 得 时 间 内 
即 可 完 水 速度 非常 决 ， 远 高 于 普通 东 发 

< 

人 
图 2 内 燕 实验 系统 。 
Fig. 2 Flash evaporation experimental Re 


@ 热电 阻 温度 计 
@ 压力 传感器 
@ 流量 计 

za 闪 


表 1 蒸发 率 的 计算 值 与 实验 值 对 比 
Table 1 Comparison between calculated and 
experimental results 
_case Tao/K pev/kPa Wa/%% mo/ 负 error/ 兄 - 


1 414.7 120 6.76 6.92 一 2.3 
2 415.6 126 6.68 6.82 一 2.0 
3 419.2 131 7.13 7.35 一 2.9 人 


3.2 内 葵 过 程 典型 曲线 Fig. 3 pr 
通过 求解 多 个 工 况 下 的 液 滴 内 燕 方 程 发 现 ， 液 
滴 温 度 、 直 径 、 质 量 、 蒸 发 速率 和 蒸发 率 在 各 工 况 下 ”3.3 不 同 初始 温度 下 的 液 滴 内 蒸 特性 
的 变化 规律 相同 、 曲 线形 式 一 致 。 以 液 滴 初始 温度 图 8~10 分 别 是 5 组 不 同 初始 温度 的 液 滴 


418 KK、 初始 直径 1 mm、 闪 蒸 压力 131 kPa 为 例 , 计 ”在 131 kPa 的 内 花 压 力 下 发 生 闪 蒸 时 ， 其 温度 、 
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图 4 液 滴 直径 变化 曲线 


Fig. 4 Evolution of droplet diameter 


mal0 tg 


图 5 液 滴 质 量变 化 曲线 


Fig. 5 Evolution of droplet mass 


mo/mg:s”! 


t/s 


图 6 液 滴 蒸发 速率 变化 曲线 


Fig. 6 Evolution of droplet evaporation rate 


燕 发 速率 和 实时 蒸发 率 随 时 间 变 化 的 曲线 . 初始 温 
度 的 不 同 仅 影响 一 开始 较 短 时 间 内 的 温度 曲线 和 温 
度 变 化 速度 ， 最 终 的 平衡 温度 是 相同 的 .初始 温度 


越 高 , 温度 降低 的 速率 越 快 . 类 似 地 , 燕 发 速率 的 差 
异 也 只 在 初始 阶段 较为 明显 ， 温 度 越 高 的 液 滴 ， 薰 
发 速率 越 快 ; 随 着 时 间 的 推移 ， 阁 发 速率 均 迅速 减 
小 趋 近 0. 由 图 10 可 以 看 出 ,初始 温度 越 高 ， 最 终 
达到 的 蒸发 率 越 大 。 在 一 定 的 内 蒸 平衡 压力 下 ， 初 
始 温度 的 提高 意味 着 过 热度 的 增 大 , 由 式 (2) 可 知 ， 
过 热度 增 大 则 传 热 系数 增 大 ， 从 液 滴 内 部 传递 的 可 
用 于 蒸发 的 能 时 增多 , 故 蒸发 率 增 大 . 图 11 则 表明 ， 
初始 温度 对 完成 内 敬 所 需 时 间 基 本 无 影响 ， 


图 7 液 滴 蕉 发 率 变化 曲线 
Fig. 7 了 Evolution of droplet flash efficiency 
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图 8 不 同 初 漫 下 的 温度 变化 曲线 
Fig. 8 Evolution of temperature of droplets at different initial 


temperatures 


3.4 不 同 内 荣 压 力 下 的 液 滴 内 蒸 特 性 

图 12~14 是 初始 温度 418 K 的 液 滴 在 6 组 不 同 
的 压力 下 进行 内 莱 时 ， 温 度 、 茹 发 速率 和 蒸发 率 随 
时 间 变 化 的 曲线 . 闪 薰 压力 越 高 ， 最 终 的 平衡 温度 
越 高 , 蒸发 速率 越 低 ,最 终 达 到 的 蒸发 率 也 越 低 . 蒸 
发 率 的 降低 是 由 于 过 热度 的 减 小 带 来 的 传 热 系数 的 
减 小 , 进而 引起 传 热量 的 减 小 所 致 但 是 由 图 15 看 
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出 燕 发 时 间 基 本 不 受 内 莹 压力 的 影响 . 
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图 9 不 同 初 漫 下 的 燕 发 速率 变化 曲线 
Fig. 9 Evolution of evaporation rate of droplets at different 


initial temperatures 
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图 10 不 同 初 温 下 的 燕 发 率 变 化 曲线 
Fig. 10 Evolution of flash efficiency of droplets at different 
initial temperatures 
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图 11 闪 蒸 时 间 与 初始 温度 的 关系 


Fig. 11 Relationship between flash duration time and initial 


temperature 
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图 12 不 同 压 力 下 的 温度 变化 曲线 


Fig. 12 Evolution of temperature of droplets at different 


evaporation pressures 
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图 14 不 同 压 力 下 的 蒸发 率 变 化 曲线 


Fig. 14 Evolution of flash efficiency of droplets at different 


evaporation pressures 
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图 17 不 同 粒 径 下 的 蒸发 速率 变化 曲线 
图 15 内 营 时 间 与 内 薰 压 力 的 关系 Fig. 17 Evolution of evaporation rate of droplets with 


Fig. 15 Relationship between Hash duration time and different diameters 


evaporation pressure 


3.5 不 同 粒 径 的 液 滴 内 蒸 特 性 
图 16~18 是 5 TS 
均 为 418 K、 闪 蒸 压力 均 为 131 kPa 的 条 件 下 菠 


线 . 可 以 看 出 ， 粒 径 的 大 小 并 不 影响 最 终 竟 用 
度 ， 只 影响 温度 变化 速度 ; 粒 径 越 坟 SN 壮 度 降低 的 
越 慢 ， 但 是 粒 径 越 大 ， 初 始 阶 司 的 族 发 速率 越 忆 
即使 如 此 , 较 大 的 液 滴 仍 需要 和 绞 长 的 时 间 完 成 燕 发 ， 
最 终 达 到 的 蒸发 率 是 相同 的 通过 图 19 可 以 看 出 ， 
闪 燕 时 间 与 粒 径 呈 线性 关系 .虽然 液 滴 的 大 小 并 不 
影响 其 最 终 的 蒸发 率 ， 但 在 工程 应 用 中 ， 考 虑 到 内 必 络 5 不 同 粒 径 下 的 花 改 率 变化 曙 线 
蒸 空间 的 大 小 有 限 ， 液 滴 停留 时 间 受 限 ， 应 选择 雾 pig de Colution of flash efficiency of droplets with different 
化 性 能 良好 的 喷嘴 ,使 液 滴 能 在 有 限 的 时 间 内 完成 ”、…、@ diameters 
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l/s 
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Fig. 16 了 Evolution of temperature of droplets with different Fig. 19 Relationship between flash duration time and droplet 
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4 结 论 


本 文 以 内 蒸 蒸 发 机 理 的 特殊 性 为 出 发 点 ， 建 立 
了 全 新 的 过 热 水 滴 闪 蒸 数学 模型 ， 并 利用 实验 结果 
验证 其 有 效 性 .计算 得 出 闪 蒸 过 程 中 液 滴 温度 、 直 
径 、 质 量 、 蒜 发 速率 和 蒸发 率 随时 间 变 化 的 曲线 , 并 
ee 
径 对 蒸发 特性 的 影响 ,得 到 以 下 主要 结 
蒸发 速率 受 由 内 向 外 的 传 热 速率 控制 . 
2) 液 滴 闪 莹 过程 中 , 温度 、 直 径 、 质 量 、 蒸 发 速 
率 和 营 发 率 在 初始 阶段 迅速 变化 ， 随 后 变化 速度 减 
组 直至 不 变 . 液 滴 闪 藻 过 程 进行 的 十 分 迅速 。 

3) 在 闪 蒸 过 程 中 , 液 滴 初始 温度 影响 其 温度 变 
化 速度 、 营 发 速率 和 蒸发 率 ， 不 影响 闪 蒸 平衡 温度 
和 闪 蒸 时 间 。. 

4) 闪 莹 压力 影响 平衡 温度 
不 影响 内 蒸 时 间 。 

5) 液 滴 粒 径 


提供 3 


、 藻 发 速率 和 蒸发 率 ， 


影响 其 温度 变化 速度 、 
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